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ABSTRACT 

We use numerical simulations of turbulent, multiphase, self-gravitating gas orbiting 
in the disks of model galaxies to study the relationships among pressure, the vertical 
distribution of gas, and the relative proportions of dense and diffuse gas. A common 
assumption is that the interstellar medium (ISM) is in vertical hydrostatic equilibrium. 
We show that the disk height and mean midplane pressure in our multiphase, turbulent 
simulations are indeed consistent with effective hydrostatic equilibrium, provided that 
the turbulent contribution to the vertical velocity dispersion and the gas self-gravity are 
included. Although vertical hydrostatic equilibrium gives a good estimate for the mean 
midplane pressure (i-')midpiane, this does not represent the pressure experienced by most 
of the ISM. Mass- weighted mean pressures {P)p are typically an order of magnitude 
higher than (P) midplane because self-gravity concentrates gas and increases the pressure 
in individual clouds without raising the ambient pressure. 

We also investig a te th e ratio Rmoi = ^H2/^hi for our hydrodynamic simulations. 



Blitz &: RosolowskyI (120061 ) showed that Rmo\ is proportional to the estimated midplane 
pressure in a number of systems. We find that for model series in which the epicyclic 
frequency k and gas surface density S vary together as k oc S, we recover the empirical 
relation. For other model series in which k and S are varied independently, the mid- 
plane pressure (or S) and i?moi are not well correlated. We conclude that the molecular 
fraction - and hence the star formation rate - of a galactic disk inherently depends on its 
rotational state, not just the local values of E and the stellar density p*. The empirical 
result -Rmoi oc (-P)midpiane implies that the three "environmental parameters" k, S, and 
are interdependent in real galaxies, presumably as a consequence of evolution: real 
galaxies trend toward states with Toomre Q parameter near unity. Finally, we note that 
i?inoi in static comparison models far exceeds both the values in our turbulent hydro- 
dynamic simulations and observed values of -Rmoh when S > IOMq pc~^, indicating 
that incorporation of turbulence is crucial to obtaining a realistic molecular fraction in 
numerical models of the ISM. 
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Introduction 



All phases of the interstellar medium (ISM) are turbulent, and this turbulence has many ef- 
fects. In the astrophysical literature, turbulence is often treated as yielding a simple addition to 
the thermal pressure, -Ptotai = pi^^ + ^turb)' "^here f^yj-b dispersion in the (one-dimensional) 

turbulent velocity, and = P/p = fk^T/fi for gas with a total number density fn and mass 
density ^n. This approach is often adopted when analyzing the stratification of interstellar gas 
clouds and the ISM as a whole, with the combined pressure gradients taken to balance the gravita- 
tional force per unit volume such that hydrostatic equilibrium is maintained by the total pressure. 
The turbulent pressure is believed to be especially important in the cold components of the ISM, 
for which observed linewidths far exceed the values of Cs inferred from excitation of atomic and 
molecular lines. 

Models of effective hydrostatic equilibrium in the vertical direction, usually assuming the 
turbulent and thermal velocity dispersions are constants independent of height, are often ap- 
plied to observatio ns of the large-scale Galac t ic ISM, and to obser v ations of the ISM in exter- 
nal galaxies (e.g. 



Lockman &: Gehman 



1991 



Malhotra 



1994 



1995 



Combes &: Becauaert 



1997 



Oiling k Merrifieldl 12001 



200e 



^; Kasparova fc Zasov 



Naravan fc .Tod l2002l : 



200 



Dalcanton et a 



For example, iNaravan &: Jog 



2004 : blitz fc Rosolowskv 



2004 



showed that the observed 



atomic and molecular disk thicknesses in the Milky Way can be fit well by assuming effective hy- 
drostatic equilibrium, and accounting for both the gas se lf - gravi ty an d the external gravitationa l 



potential of stars and dark matter. iBlitz &: Rosolowskvl (|2004l ) and iBlitz &: Rosolowskvl (|2006l ) 



(hereafter BR06) used a simplified approach to hydrostatic equilibrium in order to estimate the 
midplane gas pressure in a sample of disk galaxies, adopting a single velocity dispersion for the 
gas, treating the gravitational potential as dominated by the stars , and assuming the stellar disk's 
scale height is independent of radius. iKasparova &: Zasovl ([2003) extended the analysis of BR06 
but instead of adopting a constant scale height for the stellar disk, they assumed that the velocity 
dispersion for the stars is consistent with a state of marginal gravitational instability (with Toomre 
parameter Q^, = 1.5) for the corresponding stellar surface density. They then assumed hydrostatic 
equilibrium for all (gaseous and stellar) components separately, and computed the self-consistent 
midplane pressure, finding differences of order 30 — 40% from the simplified BR06 approach. Al- 
though widely adopted, the effective hydrostatic equilibrium model for the large-scale ISM has not, 
to our knowledge, been explicitly verified using actual turbulent flows. One of the goals of this 
work is to test this formulation systematically, using the solutions of time-dependent numerical 
hydrodynamic simulations of turbulent, multiphase gas. 

In addition to providing support against gravity, pressure also affects the phase balance in 
the ISM. For a static system at a given mean density n, changing the pressure alters the pro- 
portions of mass divided between dense clouds and diffuse intercloud medium; e.g. for cold and 
warm components in pressure equilibrium, the mass ratio of cold to warm gas is Mcoid/M„arm = 
[n/n„arm " !]/[! ~ n/ricoid] = ["-^TVarm " P]/[P ~ ^^^'coid]- The mean density itself, however, de- 
pends on pressure through the condition of vertical hydrostatic equilibrium. Turbulent pressure, as 
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it affects the response to external and self-gravity, can be expected to change both the mean density 
and the mass fractions of dense and diffuse gas. Here, we investigate these effects quantitatively. 

The fraction of ISM mass in dense gas is important from the point of view of galactic evolution, 
since this component is the immediate precursor to star formation. A recent observational study of 
external disk galaxies by BR06 identified a linear relationship between the mean ratio of molecular- 
to-atomic mass, i?moh and an estimate for the total midplane pressure oc y^S, where is the 
stellar volume density and S is the total gaseous surface density. BR06 propose that the molecular 
fractions in widely-varying types of galaxies - and hence their respective star formation efficiencies 
— are therefore determined essentially by a single parameter, the midplane pressure. To investigate 
this proposal, we use multiphase turbulence simulations in which we independently vary the input 
galactic "environmental" parameters. The observational study of BR06 focused on the dependence 
of i?moi on and S, but another important - and independent - environmental parameter is the 
angular rotation rate O (and the associated epicyclic frequency = R~'^d{Q,'^R'^)/dR). Using our 
data sets from turbulent simulations, we compare the pressure estimate of BR06 to the true value 
of the pressure, and also test how -Rmoi relates to the mean pressure measured in two different ways. 

We note that a number of recent numerical studies have investigated the formation of ISM 
structures with internal densities reaching those similar to Giant Molecular Clouds (GMCs) . Some 



studies fe.g kovama fc InutsukalbonilHeitsch et al. I bnod. bonfil : iHennebelle fc Audit 1120071 : iHennebelle et al 



2007 



2008 



Vazquez-Semadeni et al. 



20071 ) have focused on how this may occur as a consequence of 



the collision of large-scale high-velocity flows that shock and cool, becoming turbulent at the same 
time. Other studies have focused on the ability of self-gravitating instabilities to induce converging 
flows over sufficiently lar ge scales that massive, high-column density struct ures similar to observed 



GMCs are created (e.g. iKim &: Ostriker 



2001 



2007 



Li et al 



2005 



2000 ): these models include 



the galactic shear and rotati on that are important on these larg e scales, and in some cases also 
include magnetic effects (e.g. iKim. Ostriker. fc Stondl2002l . |2003| ). As spiral arms are observed to 
be strongly associated with high molecular fractions and star formation, some studies have focused 
on t he interaction between large-scale spira l shocks and self-gr avity in inducing CMC formation 
(e.g. iKim &: Ostrikeii l2002l : iKim &: Ostrikeii 120061 : iDobbd 120081 ). The details of conversion from 
diffuse to dense gas by cooling do wnstream from spiral shock fronts ha s also recently been studied 
in the absence of self-gravity (e.g. iDobbs et al. 1120081 : iKim et al.ll2008l ). Taken together, these and 
other recent studies have shown that significant quantities of dense gas form naturally as a result 
of large-scale ISM dynamical processes. Of course, dense gas in the ISM is also returned to the 
diffuse phases by the energetic inputs from star formation. In the present work, by incorporating 
feedback, we are able to evolve our models until a quasi-steady state is reached. This enables an 
analysis of the correlations among statistical properties of the system, in terms of their influence 
on the fraction of dense gas when the system has reach a quasi-steady state of cloud formation and 
destruction. 

This paper is organized as follows: In §2 we briefly summarize our numerical methods. The 
specification of model parameters and the results of statistical analysis in comparison to the vertical- 
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equilibrium approximation are presented in §3. In §4, we discuss the molecular fraction and in- 
vestigate how it relates to the ISM pressure in our models. We summerize our results and discuss 
implications for ISM structure and evolution in §5. 



2. Numerical Methods 

The analysis in this paper are based on time-dependent numerical hydrodynamic simulations 
of turbulent, multiphase, interstellar gas. Details of our numerical methods are presented in a 
companion paper (Koyama & Ostriker 2008, hereafter Paper I); here, we briefly summarize the 
model properties and parameterizations. The models we use are two-dimensional, representing 
slices through the ISM in radial-vertical {R — z) planes. We include sheared galactic rotation, 
a radial gravitational force (the centrifugal force and gravity balance in the unperturbed state, 
which assumes a rotation curve Vc = const), and Coriolis forces in the equations of motion, as 
well as gaseous self-gravity and vertical gravity representing the potential of the stellar disk. The 
gas is treated as a single fluid in chemical equilibrium, and we include (volumetric) radiative 
heating and cooling processes as a function of density and temperature appropriate to the range 
10 < T < 10"^ K. The thermal processes we incorporate include photoelectric heating from small 
grains and poly cyclic aromatic hydrocarbons, heating and ionization by cosmic rays and X-rays, 
heating by II2 formation and destruction, atomic line cooling from Hydrogen Lyman a, CII, 01, Fe 
II, and Si II, rovibrational line cooling from H2 and CO, and atomic and molecular collisions with 
grains. We adopt shearing-periodic boundary conditions in the radial direction. 

To drive turbulence, we also include a model of stellar feedback: within "HII regions" (which 
are defined by contours of the perturbed gravitational potential surrounding regions where the 
density has exceeded a specified threshold), the gas heating rate is increased by a factor 1,000. 
As a consequence, gas within these "HII regions" heats to temperatures ~ lO^K , irrespective of 
density. The detailed recipe for the feedback phenomenon is described in lPaper^ . Our aim is not 
to represent star formation feedback in a fully realistic manner, but to drive turbulence in a way 
similar to that which occurs within the dense ISM. In this sense, our feedback prescription is similar 
in spirit to sim ulations of giant rnolecular clouds i n which turbulence is applied via arbitrary forcing 



functions (e.g. IStone et al.lll998l : iMac Lowlll999l : iKlessen et al.ll200d ). Thus, our results should be 



taken as demonstrating the physical importance of turbulence to setting properties such as the 
vertical thickness of the disk, not as giving quantitative predictions for what the value of the disk 
thickness, etc., should be. 



3. Model Series and Results 



In our local disk models, three free parameters are needed to characterize the "galactic en- 
vironment": the total surface density of the gas S, the local epicyclic frequency At, and the local 
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stellar density p*. As we assume a flat rotation curve, k = \/2n where Q is the angular rotation 
rate at the center of our domain. The stellar density is used in order to specify the vertical gravity 
= -AttGp^zz. 



Following iPaper h . we study four Series of models to explore the parameter dependence of our 



results. For each Series, we hold two quantities flxed and vary a third quantity, as follows: 

• Series Q: k/T, and y/p^/T, are constant while T, varies; 

• Series K: k and y^/S are constant while S varies; 

• Series R: n/T, and p^: are constant while S varies; 

• Series S: S and are constant while k (and varies. 



Since Toomre's parameter is proportional to k/T,, Series Q and R would have constant gaseous 
Q = kCs/{ttGT,) if the sound speed Cg were constant. The Q and R series correspond to values 
of Q = 2.1(cs/7 km s~^). Assuming a constant stellar velocity dispersion, oc ^/fh, so that the 
stellar Toomre parameter (hereafter Q,,) would also have the same value for all members of Series 
Q. In all members of the R and S Series and in the E = 15.0 Mq pc~'^ models of the Q and K 
series, we take p^, = 0.14 Mq pc^'^. In the K Series, we use k = 62.4 km s~^kpc~^, while in the S 
Series we use S = 15.0 Mq pc 

This paper focuses on how turbulence affects the vertical structure of the galactic ISM. An 
important aspect of our studies is to understand how the results differ from the situation in which 
turbulence is absent. Thus, as baselines for comparison, we have two vertical non-turbulent model 
Series: one in which the gas and stellar surface densities are proportional (Series HSP), and one in 
which the stellar surface density is constant (Series HSC). These correspond to dynamical Series 
Q and K (for HSP) and Series R (for HSC), respectively. These models are one-dimensional in the 
vertical (z) direction; each model represents the asymptotic hydrostatic equilibrium state which 
develops in the absence of any stellar feedback. 

Figure [1] shows a snapshot of the gas pressure in a dynamical model from Series Q, compared 
to the hydrostatic model from Series HSP. The density and temperature are shown for the same 
snapshot in Figure 1 of Paper I. In the dynamical model, the pressure overall increases toward 
the midplane, but there are significant variations associated with structure in the gas; for the 
particular snapshot shown, there is also a high-pressure region near the left of the figure, which is 
associated with a locally-heated star formation feedback region. The hydrostatic model shows a 
secular increase in pressure towards the midplane. 
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t = 329.5 Myr 




200 400 600 aOO 1000 100 200 



X(pc) 

Fig. 1. — Left: A snapshot of gas pressure (logarithmic color scale) from Model Qll simulation. 
For comparison, the right panel shows the pressure in the hydrostatic model (HSPll) that has the 
same total gas surface density S and stellar density p^, as Model Qll. 
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Fig. 2. — Mean vertical velocity dispersion, weighted by mass. Both the thermal Cg (circles) and 
the total (thermal + turbulent) = y/^+v^ (triangles) dispersions are shown for all Series. 
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3.1. Vertical Scale Height 

We begin by examining the vertical velocity dispersion of gas in all of the model Series. Fig- 
ure [2] shows space- and time-averages (weighted by mass) of both the thermal velocity dispersion 
Cs = (Pj p)^!"^ (circles) and the combined thermal + turbulent velocity dispersion cr^ = y/c^ + v'^ 
(triangles). The four panels correspond to the Series Q, K, R, and S. In Series Q and R, the mean 
thermal velocity dispersion decreases with increasing surface density. The reason f or this is that the 



mass fraction of cold, dense gas increases with E in all of these models (see lPaper 11 ). This is because 
gravity is lower, and gas is less compressed (both vertically, and horizontally by self-gravity), at low 
S. In Series K, on the other hand, the mean thermal speed has a local minimum at intermediate S. 
Again, this can be un derstood in terms of the mass fraction of warm gas, which is largest at low and 



high S (see lPaper 11 ) in this Series; at high S, the model is extremely active in terms of feedback 
because (with constant k) the disk is quite unstable gravitationally. For all the series in which S 
is the variable parameter (i.e. Q, K, and R), the turbulent part of the total velocity dispersion 
increases with S; this is because the higher-S models have higher feedback rates, and therefore 
increasing (or flat) turbulence levels. For Series S (with constant S), the turbulence decrease as 0, 
incre ases, as h igh k stabilizes the disk and prevents gravitational collapse and feedback (see Figure 



11 in 



Paper J ). For all series, the (mass- weighted) turbulent vertical velocity dispersion approaches 
or exceeds the (mass-weighted) thermal velocity dispersion for some part of parameter space, so 
that turbulent support of gas in the vertical gravitational field is expected to be important. 

Next, we measure (for all Series) the vertical scale height, using the following averaging: 



g,ve = W ^" (1) 
y 2-^ all zones P 

where z is the vertical coordinate relative to the midplane. We further average the values of i^ave 
over time. In order to test whether the velocity dispersion can be used to obtain an accurate 
measure of the scale height, we also compute "estimated" vertical scale heights defined as: 

^ ^ J: ^1 1 

^/2^GS + [(GS)2 + 2G/5,a2]V2 A + [^2 + 1] ^' 

' (S - 0) 



(3) 

(p*-0). 



V87rGS 

This formula (see Appendix for derivation) accounts for both gas self-gravity and stellar gravity; 
the limiting forms are for negligible gaseous and stellar gravity, respectively. 

In equation Q, ^ is a dimensionless factor that measures the relative densities of the gaseous 
and stellar disks, 
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Fig. 3. — Disk scale heights, for all hydrodynamic and hydrostatic models. Open boxes denote 
the directly- measured scale height (see eq. [1]) for all hydrodynamic Series. Filled boxes show the 
measured scale height for corresponding hydrostatic models (HSP for Series K and Q, HSC for 
Series R). Open circles and triangles show the estimated scale heights (see eq. [2]) using thermal 
and thermal plus turbulent velocity for az, respectively. The bottom part of each panel shows the 
ratio of estimated scale heights to direct measurements. 
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The latter expression treats the stellar disk as an isothermal self-gravitating equilibrium, with = 
c1 ^/(ttGS*), and shows that A ~ Q^jQ (assuming that vertical and radial velocity dispersions are 
proportional). The formula ([2]) may be thought of as an extension of the usual non-self-gravitating 
scale height formula to account for the gravity of the gas. Since A > 0, the correction factor 
depending on A is always < 1. If the gas disk is much more gravitationally unstable than the 
stellar disk (A ~ Q^jQ ^ 1), the correction factor is large; otherwise the correction factor is 
order-unity. 

Figure [3] shows the measured (i^ave) and "predicted" (-ffest) disk scale heights for all series of 
hydrodynamic models. For -ffest; we show results using for Oz either the thermal velocity dispersion 
{(Jz = Cs] subscript c) or the total velocity dispersion (cr^ = + v^; subscript c + Vz)- To show 
how turbulence contributes to setting the disk thickness, -ffave is also shown for the hydrostatic 
models. The difference between H^ve in hydrostatic and hydrodynamic models can be quite large, 
up to a factor 10 in some cases. We note that .ffave of the hydrostatic models {filled boxes) differs 
from Hcst,c {open circles) because the mass- weighted mean sound speed differs for hydrostatic and 
hydrodynamic models. 

Overall, Figure [3] shows that the estimate for scale height Hcst,c+vz that includes turbulence 
traces the measured -ffave quite well, for all the Series. The difference between -ffest,c and -ffest,c+vz 
increases with increasing E, with quite large differences for some of the models in Series Q and R. 
This indicates that high surface density disks are supported largely by turbulent velocities, in these 
cases. To facilitate comparisons between estimated and measured value of the scale height, in the 
lower part of each panel we also show the ratios -ffest,c/-f^ave {circles) and -ffest,c+i)z/-f^ave {triangles). 
At low values of S in Series Q, K, and R, both estimates of H exceed the true measured value. 
It is notable that where the turbulent contributions are large, at high S in Series Q, K, and R, 
the estimated and measured disk thicknesses are in quite good agreement (within ~ 10 — 20%). 
Thus, we conclude that if measurements of the vertical velocity dispersion together with the gaseous 
surface density and stellar surface density can be made observationally, they can be combined to 
yield an accurate estimate of the gas disk's thickness. 



3.2. Gas Pressure 

The gaseous pressure, like the scale height, is often difficult to measure directly. As a conse- 
quence, other proxies are often used to obtain an estimate of the value of the pressure, with an 
assumption that vertical equilibrium is satisfied. Here, we test how well such pressure estimates 
agree with the directly-measured pressure, for our multiphase turbulent models. 

Figured] shows for all models in all hydrodynamic Series the average gas pressure. We consider 
two different ways of averaging: weighting by mass {P)p {open box), and weighting by volume 
(-P)midpiane {opcn circlc). The value {P)p is interesting because it characterizes the value of pressure 
experienced by the average atom or molecule, whereas (P)midpiane is interesting because it represents 
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the pressure in the diffuse (non-self-gravitating) part of the ISM that is closest to star-forming 
regions. 

The mass- and volume- weighted averages are defined by the following: 

(-P)midplanc = / ^ J^- (6) 

For {P)p, all zones in the domain are included, while for (-P)midpiane) the subscripts and ^ + 1 
indicate that only zones in the two horizontal planes closest to the midplane are included. Time 
averaging is applied in all models after the above space averaging. We also show the same pressure 
averages for the hydrostatic Series {filled box and filled circle). Interestingly, in the hydrodynamic 
models {P)p always exceeds (-P)midpiane by a large factor ~ 10. This indicates that self-gravity is 
important in increasing the pressure above the "ambient" value, for much of the gas. Pressures 
cannot exceed the ambient midplane value without horizontal gradients, which are balanced by 
the gravity within individual clouds (see Fig. [1]). In the hydrostatic models, {P)p {filled boxes) 
is generally slightly below (P) midplane {filled circles), because the pressure at the midplane is the 
maximum within any system, and weighting by mass includes lower-pressure gas which reduces the 
average. (Note that for the hydrostatic models, there are no horizontal gradients in any quantities; 
see Figure [TJ) Except in the most active disks, the mass- weighted averages for the hydrostatic 
models are close to the midplane values for the hydrodynamic models. In Figure H] we also display 
the pressure estimate of BR06 {solid line) defined as: 



Pbr = ^vy^2G^,, (7) 

where v = 8 km/s is adopted. This line falls between {P)p and (-P)midplano for all the hydrodynamic 
Series. 

For hydrostatic Series HSP (shown in the Series Q and K panels), the slope of the midplane 
pressure is close to that predicted by equation d?]), while being offset to lower P by a factor 2-3. 
The difference in slope is because the medium has multiple phases, rather than a single phase at a 
given thermal sound speed. The offset is because (i) much of the mass in the hydrostatic models 
is at low temperatures, for which the sound speed is well below 8 km/s, and (ii) equation d?]) does 
not include the gaseous vertical gravity, which is comparable to the stellar gravity when vertical 
velocity dispersion is low and the disk is very thin (see below). These effects push P in opposite 
directions, and hence partially compensate each other. For hydrostatic series HSC (shown in the 
Series R panel) , the prediction of equation ([7]) departs significantly from the slope of the midplane 
pressure results, because in the HSC series (which has p* constant) vertical gravity is strongly 
dominated by gas rather than the stellar component at large S. 

In ^3.H we defined an "average" vertical equilibrium using the total surface density and the 
total vertical velocity dispersion, and showed that this could yield an accurate measurement of 




Fig. 4. — Gas pressure averages for all Series. Open boxes show mass-weighted averages and open 
circles show the midplane pressure, for hydrodynamic models. Filled boxes and filled circles show 
the same for hydrostatic models. The pressure estimate of BR06 is also indicated solid line. 




Fig. 5. — Measured, fitted, and estimated gas pressures. Open boxes and circles are the same as in 
FigureHl Solid and dashed lines are the corresponding fits to equation (jlOp . Our vertical-equihbrium 
midplane thermal pressure estimate (eq. [9]) is plotted as filled triangles. 
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the disk thickness. The same model (see Appendix) can be used to estimate a midplane "effective 
hydrostatic pressure," which we can compare to measured values. If H is the scale height, then in 
equilibrium the mean midplane gas density is pQ = T,/{^/2TrH). Using the total velocity dispersion, 
the predicted total gas pressure at the midplane is then given by Po.tot = c^/^O; which using equation 
([2]) gives 

Po,tot = = EfGS+r(GS)2 + 2Gp,a21^/' 



The expression ([8]) corresponds to an extension of the pressure estimate formula of BR06 using the 
inverse of the ^-dependent correction factor that appears in the scale height estimate ([2]). This 
correction factor is > 1. 

Equation ([8]) gives an estimate of the total midplane pressure, but the thermal pressure should 
represent only a fraction (c^) / {c^ + v^) = 1 — {v1) /a1 of Po.tot, where (w^)^^^ is the mass- weighted 
RMS turbulent velocity dispersion in the vertical direction. Thus, our estimate of the mean thermal 
pressure at the midplane is 

,2\ 



Po,th = {cl)^=^ = S^y2G^(yl + (9) 

In Figure O we compare the pressure estimate from equation ([9]) {filled triangles) with the measure- 
ments of midplane pressure. At large E, the agreement is quite good, while at low S the estimated 
midplane pressures lie slightly above the measured values. This behavior is similar to our results 
for estimated scale heights, which were in best agreement with the measured -ffave at large S (where 
the dense gas dominates the mass, and the velocity dispersion is turbulence-dominated). 

For all the Series in which S is the independent variable, we have fit the measured gas pressure 
to the formula: 



P/kB = Dj ^^ ( ^„-2 ) • (10) 



p* / ^ 

M© pc-3 VMq PC- 

For {P)p and (P) 

midplane) we find the respective coefficients are Dp — 1.3 x 10^ K cm and 
Djnidpianc = 1.1 X 10'^ K cm"^, respectively. The largest and smallest surface density models are 
excluded in the fits. The results of the fits are displayed as solid and dashed lines, respectively, in 
Figure [5l To compare with the BR06 formula, we also fit {P)p and (-P)midpiane to P = CT,v\J2Gp^ 
with w = 8 km/s. We find Cp = 3.6 and Gmidpiane = 0-3; i.e. the BR06 formula for pressure yields 
values that are typically a factor ~ 3 larger than our measured midplane pressures, and a factor 
~ 4 below the mass-weighted mean values of pressure. As noted above, the mass-weighted average 
pressures are about ten times larger than the midplane pressures; this is evident in the ratio of the 
fitting coefficients. 

Finally, we note that for most models (except at low S), the measured midplane pressure 
exceeds the maximum pressure of the warm neutral medium, Pu,,max/^B = 5.5 x 10^ K cm^^ for 
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our adopted heating and cooling functions. Dense clouds that are externally confined by the warm 
medium cannot have pressure exceeding P^^max unless they are internally stratified (implying they 
are self-gravitating); thus, -Pio,max is the largest the midplane pressure could be in the absence of 
self-gravity. Equation ([9]) can be solved for S in terms of the midplane value of -Po,th- The maximum 
surface density for an atomic-only disk without self-gravitating clouds is then obtained by setting 
-Po.th Pw,ma.x, with the result S (Pw, ma.x / G)^^"^ az / Cg times a function of A that varies between 
0.3 and 0.6 for A = 0.1 — 1. Assuming az/cg ~ \/2 and taking P«,,max/^B = 5.5 x 10'^ K cm~^, 
the maximum surface density for a pure-ator nic disk is ~ l O Mf?) 



PC ^; this is consistent with the 



saturation levels for HI gas observed e.g. bv lWong &: Blitzj (|2002l ). Since the measured midplane 



pressure is a volume- weighted sum of the pressures in different phases, a mean value exceeding 
-fw,max implies that self-gravitating dense clouds occupy a non-negligible fraction of the midplane 

volume, fv = {Mdense/MdiSusc){PdiSusc/ Pdcnsc), with (P)midplano = (-Pdense " -^diffuse ) /v + -PdiflFuso • In 

the next section, we turn to a discussion of the relationship between the dense-to-diffuse mass ratio 
and global parameters. 



4. An Application: Molecular Mass/Pressure Relations 

In this section, we explore relationships between the dense gas fraction and "environmental" 
conditions, including the gas pressure and the gas surface density. We are motivated by observations 
that show high molecular fractions in environments - including spiral arms and galactic center 
regions - where both the total gas surface density and stellar density are high. In particular, BR06 
found for a number of disk systems that the mean ratio of molecular-to-atomic mass scales nearly 
linearly with the pressure estimate Pbr defined in equation ([7]). Although our turbulent, multiphase 
simulations show that Pbr in fact overestimates the pressure of the typical volume element and 
underestimates the pressure of the typical mass element, Pbr nevertheless systematically increases 
in a similar way to both (P)midpiane and (P)p. Thus, it is interesting to test how the dense-to-diffuse 
gas mass ratio depends on the true values of pressure. In addition to empirical results suggesting a 
relation between mass ratio and pressure, there are theoretical reasons that the mass ratio should 
depend on the mean gaseous surface density. For example, if atomic gas is converted to molecular 
clouds through gravitational instabilities on a timescale tfomi ~ f'"Hi/(GS), and molecular clouds 
are destroyed by star formation on a timescale tdest j equating cloud formation and destruction rates 
implies Mn^/Mui = tdest/^fomi) which is oc E if the HI velocity dispersion and cloud destruction time 
are relatively constant. Thus, it is interesting to explore dependence of Mu^/Mm on the surface 
density - which appears in both the effective hydrostatic pressure and the rate of self-gravitating 
instabilities. 
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4.1. Molecular Gas 



Although our numerical model does not directly include formation/dissociation processes of 
H2, we can nevertheless relate our results to observed gas phases in an approximate way, using 
density as a proxy. Namely, we expect gravitationally bound dense clouds at n > 100 cm^^ to 
consist primarily of H2, whereas diffuse gas at lower densities consists primarily of HI. We argue 
for this approximate identification based on the formation/dissociation equilibrium condition for 
H2 molecules, which includes photodissociation and cosmic ray dissociation, and formation on dust 
grains: 



pump 



+ 



CcR)^/mol 



(11) 



(jTielens fc HollenbachI Il985l ). Here, /moi = 2n(H2)/n is the molecular fraction, and 1 — /moi = 
n(HI)/n is the atomic fraction. The FUV dissociation rate is limited by shielding, which depends 
on the optical depth in H2 lines and the extinction. Formation on grains depends on the sticking 
probability. The details of the terms involved are listed in Table [H We adopt FUV field strength 
Go = 1.7, gas and dust temperature T = 10 K, and cosmic-ray ionization rate of hydrogen atoms 
= 1-8 X 10~^^s~^. For any total hydrogen column A^h and volume density n, we can solve to 
obtain /moh the molecular fraction. Figure [6] shows, in the n — N plane, the boundary {solid line) 
between the predominantly-atomic and predominantly-molecular regimes, which we define by the 
locus of points for which /moi = 0.5. 

At any density, we can also define the Jeans length Lj = Cs{tt/Gp)^^'^ , where = fceT/^ 
(we adopt T = lOK). This defines a corresponding total column of gas, Nn = nLj, that could 
be expected to be gravitationally bound. The boundary between gravitationally unbound (low n 
and A^) and bound (high n and N) gas, based on this criterion, is shown in Figure [6] as a dashed 
line with n oc A^^. Note that if instead of Lj we had chosen as a l e ngth scale t he diameter D 
of sphere containing mass equal to the Bonnor-Ebert (jBonnoii Il956l : lEbertl 119571 ) critical mass, 
Mbe = 1.182c^/(G^p)i/2^ ^j^gn ^ 0.74Lj. This would shift the unbound/bound line in the 
log(A^) - log(n) plane to the left by log(0.74) = -0.13. 

We note that the gravitational binding criterion discussed above considers only support by 
thermal pressure. Turbulence can lend further support against gravity, and this is particularly im- 
portant for molecular gas, which is quite cold. For example, if we considered turbulence-supported 
clouds with velocity dispersion following the observed linewidth-size relation of Galactic GMCs 
(jSolomon et al.lll987l ). then the column density separating gravitationally bound from unbound 
regions would have a constant value equal to half of the mean observed CMC column, amounting 
to A^H = 7.5 X lO'^^ cm A higher nor malization for the linewidth-size relation (as occurs in galac- 



tic center regions; see lOka et al.l[200ll l would further shift the unbound/bound limit to larger A^. 



Thus, moderate-density molec ular gas can in pr inciple be gravitationally unbound under conditions 
of sufficiently high turbulence (lElmegreenI 1 1 9931 ) . For our current simulations, however, turbulence 
levels are not this high (see discussion below). 
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Bounded HI 



Fig. 6. — Phases in the density-column density plane. The solid line divides the area into 
predominantly-HI at low nn and TVh, and predominantly-H2 at high nn and Nn, adopting molecule 
formation and destruction processes as described in the text. Shown as a dashed line is N = nLj, 
where Lj is the Jeans length at n for T = lOK gas. The left and right hand sides of this line are 
gravitationally unbound and bound, respectively. Dotted hues show the typical resolution limit of 
the simulations Ax = 2 pc, and a maximum cloud scale of 100 pc. 
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We have denoted the three different regions in the log(A^) — log(n) plane according to their 
expected chemical and gravitational properties. The crossing point of the two separation loci is at 
n ~ 125 cm~^ and iVn ~ 1.2 x 10^^ cm ~^ (Av ~ 0.75 mag), with corresponding local Jeans length 
of Lj = 3.1 pc. This size is in fact slightly larger than the typical resolution limit of our simulations, 
Ax = 2pc; we show this limit in Figure [6] as a dotted line, with regions to the right resolved and 
those to the left below the resolution limit. The resolution limit crosses the HI/H2 separation 
curve at n = 180 cm^'^. Because the resolution limit falls at larger than the bound/unbound 
separation nearly everywhere in the molecular domain, all zones at a given density that are resolved 
and molecular would also be gravitationally bound. In practice, clouds do not exceed ~ 100 pc in 
crossection; we have marked this limit in the Figure as a dotted line. 

According to the limits shown in Figure El any resolved regions in our simulations at n > 
100 cm ~^ would be molecular. This is a conservative definition, since it omits some gas between 
n ~ 10 — 100 with Ah > 10^^ cm that could be molecular. However, gas at these densities 
could also be in the cold atomic phase (which extends down to ncoid,min = 8.6 cm for the 
cooling curve we adopt); we choos e the str icter definition. We note that when the virial ratio 



('--^kinetic/gravitational energy; see iPaper J ) is measured for gas in the range n = 1 — 100 cm 
(most of which is at 10cm~^ < n), the values are well above unity - implying that gas parcels in 
this density range are mostly found in non-self-gravitating regions with low surrounding column 
densities, to the left of the unbound/bound curve. This suggests that in practice very little high- 
column density gas that would be molecular is missed when we set the minimum threshold at 
100 cm From the point of view of dynamics, this is because the density rises whenever any region 
becomes gravitationally bound, so low-density regions at high column are rapidly depopulated. We 
also note that the formation time in dense, cold regions is expected to be short (~ 10^ years 



from iGlover &: Mac Low I l2007bl ). because supersonic shocks increase the density above ambient 



values and accelerate the molecule formation process, which occurs at a rate nRf. 



4.2. Molecular Mass-Pressure Relation 

Following the discussion in the previous section, we adopt a working definition of molecular 
gas as that at tih > 100 cm Atomic gas therefore consists of the lower-density complement, 
including both what would be observable as warm and as cold HI in 21 cm emission. The mass 
ratio of molecular to atomic hydrogen is then defined as 

_ M{n> lOOcm-3) 
^'""^ = M(n< lOOcm-3)' ^ ' 

where we apply space- and time-averages before taking the ratio. 

Figure [7] shows Rmo\ as a function of Pbr = '^v^/2Gpl for all hydrodynamic and hydrostatic 
Series (we use v = 8kms~^ as in BR06). We also show as a solid line the empirical fitting formula 




Fig. 7. — Mean molecular-to-atomic mass ratio i?moi as a function of Pbr (see eq. [7|), shown as 
open boxes for all hydrodynamic models. Filled boxes show the results from hydrostatic models. 
The solid line shows the empirical fit from BR06 (eq. [T3]) . 
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from the observational study of BR06 (see their eq.l3): 



mol 



nu II l0.92 



4.3 X 104 K cm -3 



(13) 



Interestingly, we find that our results for i?moi follow the empirical result for some but not all series. 
In particular, the models in Series Q and R - which have values of that scale with S in such a way 
as to keep the gaseous Toomre parameter constant - are close to the BR06 fit. The models in Series 
K, which have constant k and therefore high (or low) values of k/S where Pbr is low (or high, 
respectively), do not follow the empirical result of BR06, but instead show a ratio i?moi that is near 
unity independent of Prr- This has two interesting implications. First, our models with oc S 
have similar behavior to real galaxies, indicating that real systems evolve (by converting their gas 
to stars) in such a way as to have Toomre parameter fall within a limited range of values. Second, 
because the K Series departs from the BR06 result, our models suggest that the molecular fraction 
does not have a one-to-one relationship to the effective pressure parameter Pbr- Comparing series 
Q and K which have the same S and p^, -Rmoi increases with increasing k. For example. Figure 8 
shows that i?moi increases by factor 2.4 when Q (and k) increases by factor 2^/2, for the highest-S 
(S = 42Mq pc ~^) model. For the S = 21Mq pc model, the Rmoi increase is 60% for an 
increase by a factor ^/2, comparing the Q and K series. Series S, which varies k at a given value of 
S and also shows departures from the empirical i?moi vs. Pbk relation. 

Given that molecular gas in our models is primarily found in gravitationally-bound systems, 
it in fact makes sense that the molecular fraction should not have a one-to-one relationship to the 
parameter PbRj since Pbr does not include any effects of galactic rotation. Galactic rotation and 
shear are crucial for regulating the large-scale gravitational instabilities that create giant molecular 
clouds in real galaxies as well as in our models, so we believe that the molecular-to-atomic ratio must 
intrinsically be sensitive to environmental factors that are not captured in Pbr- Thus, if observed 
galaxies do show a one-to-one relation between -Rmoi and -PbRj it implies that the environmental 
parameters k, S, and p^, are not all independent in real systems. 

In Figure [HI we show i?moi as a function of the surface density for all of our model Series. The 
behavior is similar to that shown in Figure [7] because Pbr depends monotonically on S for all our 
Series: Pbr 5]^ for Series Q and K (which have p^^ oc S'^), while Pbr 5] for Series R (which has 
yO=K = const.). From both Figures Figure [7] and [8l it is evident that the hydrostatic models {filled 
boxes) generally have a much larger molecular component than both the hydrodynamic models 
and the empirical results, except at low gaseous surface density. This indicates that turbulence is 
essential for determining the phase balance between diffuse and dense gas in the ISM as a whole. If 
the ISM were a static system, it would be overwhelmingly molecular even at fairly moderate values 
of S and p*. In real galaxies, turbulence limits gaseous settling into the midplane and the extreme 
self-compression that would otherwise ensue. 

While we have argued that the molecular content of a galaxy cannot (in general) be predicted 
solely from S and p* because self-gravitating horizontal contraction is also responsive to the local 
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Fig. 8. — Molecular-to- atomic mass ratio -Rmoi vs. total gaseous surface density S. Open boxes 
show the results from hydrodynamic models, and filled boxes show the hydrostatic model results 
(HSP for Series K and Q, HSC for Series R). The solid line indicates the empirical result from BR 
(eq. [T3l using eq. [7|). 




Fig. 9. — Molecular-to-atomic mass ratio -Rmoi vs. measured mass-weighted mean pressure {P)p 
and midplane pressure (-P)midpiane foi' all tiydrodynamic models. Lines show linear fits (see text). 
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rotation and shear rates, it still is plausible that the molecular fraction should reflect the true 
mean pressure in the ISM. If molecular gas is collected in self-gravitating clouds, then since their 
internal pressure is higher than ambient levels, an increase in the molecular fraction should go hand- 
in-hand with a higher mass-weighted mean pressure {P)p. At the same time, ambient midplane 
pressures (P)midpiane increase when the total gas surface density increases, and (provided that k is 
low enough) larger S also renders the disk susceptible to gravitational instabilities that would form 
dense, bound clouds and increase the molecular fraction. 

We explore these ideas by plotting in Figure [9] the molecular-to-atomic ratio against our two 
measures of mean gas pressure, {P)p and (-P)midpiane- We also fit the combined results for all 
Series to single linear relations. These fits, overplotted in Figure [9] as solid and dotted lines, are 
i?moi = {P)p/[7-G X lO^cm-^KfeB] and R^oi = (-P) midplane/ [6. 7 X lO^cm-^K/sB]. For Series Q 
and R, the fits using {P)p are quite good, and the fits using (-P)midpiane are also fairly good (Series 
Q and R also show better agreement with empirical results than the other Series). For Series K, 
the fit using {P)p is reasonably close to the models results, but the fit based on (P) midplane fails in 
a similar manner to that shown in Figure ([7]) and discussed above. The basic reason for this is that 
the midplane pressure, either directly measured or estimated using equation Q, increases with 
increasing S. However, the molecular-to-atomic ratio for Series K does not strongly and secularly 
increase with S due to the differences in rotational effects in this constant-K Series compared to 
the other Series. At high values of T, in Series K, the disk is extremely gravitationally unstable 
overall, and as a consequence is more active in producing feedback than other models at the same 
S. As a consequence, a smaller fraction of the gas mass ends up being in the dense phase than in 
Series Q and R. Overall, we conclude that .Rmoi is indeed well correlated with the mass-weighted 
mean pressure, {P)p, as (almost definitionally) is expected. The measured mean midplane pressure, 
which is more closely related to simple vertical-equilibrium pressure estimates, is less well correlated 
with i?moi when environmental parameters k, and S are all independent. 



5. Summary and Discussion 

We have used numerical simulations of turbulent, multiphase, self-gravitating gas orbiting in 
the disks of model galaxies to study the relationships among pressure, the vertical distribution of 
gas, and the relative proportions of dense and diffuse gas. In particular, we compare the results on 
vertical stratification obtained from space-time averages of fully-dynamic - and often turbulence- 
dominated - systems with simple estimates based on single-component effective hydrostatic equilib- 
ria. We also investigate how vertical-equilibrium estimates for the pressure compare with measured 
mean values of the pressure in our models. Empirical studies by BR06 have identified a linear 
relation between the molecular-to-atomic mass ratio i?moi and a midplane ISM pressure estimate, 
Pbr ^y/p*- We study the origin and implications of this relation by testing the correlations 
among iimoh -PbRi and the directly-measured midplane and mean pressures in our models. 

Our chief conclusions, and their implications, are as follows: 
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1. The average disk scale height is weh represented by estimates that assume hydrostatic 
equihbrium and an effective total pressure based on the total (thermal + turbulent) vertical velocity 
dispersion (see Fig. [3] and eq. [2]). Thus, provided that gas surface densities, vertical velocity 
dispersions, and stellar density can be measured, an accurate estimate for the disk thickness can 
be obtained. 



Hydrostatic equilibrium with an effective turbul ent pressure is commonly assumed iii both 



Galactic and extragalactic observational studies (e.g. 


Lockman &: Gehman 


1991; 


Malhotra 


1994. 


1995; 


Combes & Becauaert 


1997; 


Olhne: & Merrifield 


2000; 


Naravan &; Joe 


2002; 


Dalcanton et al. 


2004; 


Blitz & Rosolowskv 


2004, 


2006; 


Kasparova & Zasov 


2008 


) , but to our knowledge the relations 



that are generally adopted have not previously been tested with direct numerical simulations. 
Our hydrodynamic studies demonstrate that for determining the scale height H, the effective 
hydrostatic equilibrium assumption is indeed sufficient, even when turbulent support far exceeds 
thermal support (and provided that magnetic effects are sub-dominant; see below). Thus, measured 
disk thicknesses in edge-on disk galaxies could in principle be used to determine the unobservable 
vertical velocity dispersion, and measured line-of-sight velocity dispersions in face-on galaxies could 
be used to determine the unobservable disk thickness. 

The basic reason the hydrostatic formula can be used to obtain an accurate measure of H is 
that what is really being equated is the total vertical momentum flux p^k^T / fi + v'^) averaged over 
the midplane, and the total vertical weight of the ISM, / dz pQz ^ pA'KG{p + p^:)H'^, averaged over 
the horizontal direction. Provided that the time-averaged value of the momentum per unit volume 
in the midplane does not change, momentum conservation including gravitational source terms 
demands that the difference between vertical momentum flux and vertical weight must be zero, 
independent of details of the dynamics. The formula ~ cr^/[47rG'(p* -|- T,/ H^/2Tr)] is therefore 
fundamentally an expression of momentum conservation. 

2. Mass-weighted mean pressures {P)p in our hydrodynamic models significantly differ from 
the mean midplane pressure (-P)midpiane5 while these quantities are quite similar to each other in our 
static comparison models. Typically, the hydrodynamic models yield values of {P)p an order of mag- 
nitude larger than (-P) midplane- The difference can be attributed to self-gravitating condensation, 
which makes concentrated clouds with high internal pressure rather than a horizontally-uniform 
gas distribution with more moderate pressure. 

Simple estimates of the pressure based on vertical hydrostatic equilibrium fall between mass- 
weighted and midplane values, with the formula used by BR06 (see our eq. [7]) comparable to the 
geometric mean Pbr ~ \/ (-P)p(-P) midplane- A single-component estimate for the midplane thermal 
pressure that accounts for self-gravity and the mean thermal and turbulent velocity dispersions 
(see eq. [9]) follows the measured midplane pressure fairly closely, especially at high S. Thus, if 
turbulent and thermal vertical velocity dispersions can be measured directly (for face-on galaxies), 
a good estimate of the midplane total or thermal pressure can be computed via equation ([8]) or 
([9]). For an edge-on system in which the scale height is measured, the midplane total or thermal 
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pressure can be estimated as S/(i?-v/27r) x (erf) or (cf ). Midplane pressure estimates based on large- 
scale observables that assume hydrostatic equilibrium can be quite accurate, but this depends on 
an accurate measure of the vertical velocity dispersion or vertical thickness. Even if the velocity 
dispersion is not known, the relative midplane pressures of different regions within a galaxy (or 
from one galaxy to another) can be obtained using the hydrostatic formulae, provided the variation 
in the (unknown) velocity dispersion within the observational sample is small compared to the 
variation in the stellar volume and gaseous surface densities. Midplane pressure estimates made in 
this way should not, however, be treated as a proxy for the pressure in the typical mass element, 
{P)p, which can be much larger than the pressure in the typical volume element. 

3. Based on calculations of molecular abundance as a function of hydrogen volume density n 
and column density N combined with the resolution and measured virial ratios in our simulations, 
we adopt a working definition of gas at n > 100 cm ~^ as "molecular" and n < 100 cm ~^ as 
"atomic". We then investigate the ratio -Rmoi = M-n^/Mm for all our models. We find that Series 
Q and R, which have rotation rate oc S, show correlations between i?moi and Pbr (or i?moi and 
S) that are similar to the empirical result reported by BR06, Rmoi -Pbr oc Sy^. On the other 
hand, Series K and S, in which S and do not vary together, depart from the empirical relation 
-Rmoi Pbr- 

We conclude that (i) the molecule fraction inherently depends on the rotational state of a 
galactic disk, not just on the local values of the stellar volume and gaseous surface densities p^: 
and S; and (ii) the empirical relation -Rmoi Per identified by BR06 implies that the third 
"environmental parameter," the epicyclic frequency k = (assuming a flat rotation curve), is 
not independent of p^: and S in real galaxies. This dependence can be accomplished by evolution: 
for example, disk galaxies may convert gas into stars until they reach a state in which the Toomre 
parameter oc k/S approaches a critical value. 

4. We have tested the correlation between i?moi and the measured pressures in our models, 
(^)midpianc and {P)p, and find a good correlation in all Series only for the latter. The correlation 
between {P)p and i?moi is potentially useful as a way to estimate the typical internal pressure within 
gravitationally-bound regions when only the total molecular-to-atomic mass is easily accessible, as 
for low-resolution observations. This internal pressure is important in the small-scale aspects of 



star formation such as determining the IMF (jMcKee &: Qstrikeij 120071 ) , as well as in molecular 



chemistry. The lack of correlation between -Rmoi and (i^) midplane in Series K implies that the 
molecular content cannot in general be predicted solely from S and p=K (i.e. without knowledge of 
k), as noted above. This reflects the fact that the formation of self-gravitating clouds is regulated 
not just by gravitational processes and pressure, but also by angular momentum. 

5. For our non-turbulent comparison models, we find that -Rmol far exceeds observed values. 
This indicates that turbulence is essential to setting the observed phase balance in the ISM. Recent 
theoretical investigations of the origin of Kennicutt-Schmidt laws have focused on the dependence of 



star formation rates on the molecular, rather than total, gas surface density (e.g. iNaravanan et al 
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20071 : iRobertson fc Kravtsovl |2008| ) . Since turbulence is crucial in determining the abundance of 
dense gas, in simulations that aim to compute this abundance realistically it is necessary to incorpo- 
rate the feedback effects that drive turbulence, and to run on a fine enough mesh (or with sufficient 
SPH particles) that the turbulence is well resolved. While technically challenging in global disk 
models, local models may offer a more immediate route to this goal. 

Caveats - The models analyzed in this paper are subject to a number of limitations, which 
could potentially affect some of our conclusions. The chief limitations of the simulations are that (i) 
they are two-dimensional, representing cuts in the R — z plane, rather than three-dimensional; (ii) 
we have adopted a very simple model to implement turbulent driving as a star formation feedback 
effect from HII regions, and we have not included other drivers of turbulence such as supernovae, 
spiral shocks, and shear instabilities; (iii) we have not included magnetic fields (or cosmic rays). 
We intend to pursue these extensions in future work. 

Inclusion of magnetic fields and altered turbulent driving would certainly affect the specific 
quantitative findings for ifave, (-P)midpiane) {P) p-, and -Rmol in our models. We believe, however, 
that the results we have emphasized regarding physical relationships are robust. In particular, 
with appropriate modifications to include magnetic stresses, the time-averaged vertical momentum 
flux through the midplane must still equal the time-averaged vertical weight if the mean vertical 
momentum is conserved. This can be used to predict the total midplane pressure (including the 
magnetic pressure) and H given the values of E, p^, and the thermal, turbulent, and Alfven 
velocities. Thus, we anticipate that inclusion of magnetic fields and alternate turbulence sources 
would not fundamentally alter the conclusion that reasonable estimates of scale heights can be 
made using observable quantities even in highly-dynamic systems. 

Further, we expect that our conclusions regarding the presence or absence of correlations 
between i?moi and {P) p or (P)midpiane would continue to hold in models that include additional 
turbulence sources and magnetic fields, although the details of correlations might change. Namely, 
angular momentum inherently must be important in permitting or preventing formation of dense, 
self-gravitating clouds. Our present models account for angular momentum effects, and show that 
Rmo\ does not in general have a one-to-one relationship with (-P)midpiane or Sy^; we expect this 
result would carry over into any model that incorporates sheared background rotation of the galactic 
disk. Thus, if a one-to-one relationship between i?moi and indeed exists empirically, it implies 

that S, p*, and k are not all independent quantities in real galaxies. 

We are grateful to the referee for a number of comments that have helped improve our presen- 
tation. Numerical computations used in this project were carried out on the OIT High Performance 
Computing Cluster, and the CTC cluster in the Department of Astronomy, at the University of 
Maryland. This work was supported by grant NNG05GG43G from NASA. 
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Table 1: Processes and parameters for H2 formation/dissociation 



Formation of H2 on dust grains 


Rf 


= 6 X W-^'^{T/300f-^S{T) cm^s-i 


[1] 


Sticking probability 


S{T) 


= [1 + 0.04(T + Tdf-^ + 2y. lO^^T + 8 x lO^^T^]-! 


[1] 


Photo dissociation rate 


-Kpump 


= 3.4 X 10-i°Go/3(r)exp(-2.5^„) s^^ 


[1] 


Self-shielding function 






[1] 


Optical depth 


T 


= 1.2 X m~^^fNn5v^^ 


[1] 


Cosmic-ray dissociation 


''CR 


= 2.29C?R 


[2] 


Turbulent line broadening 


Svd 


[3] 


Visual attenuation 




= iVH/1-5 X lO^i cm-2 





[1] iTielens fc HoUenbachl (|l985h 
[2] iHoUenbach fc McKeel (|l989h 
[3] ISolomon et all (|l987l ) 



- 30 - 



A. Vertical Equilibrium with Stellar and Gas Gravity 

The vertical momentum equation, when averaged over a horizontal plane, is given by 



dV"^ dz\ ' Stt 47r/ \ dz 



(see e.g. iPiontek &: Qstriked 120071 ) . Here, B is the magnetic field and <I> is the total (stellar 
plus gaseous) gravitational potential. In steady state {pv^) is time-independent, so if we neglect 
magnetic fields and assume that p, f^, = P/ p, and d^/dz are statistically independent quantities, 
we obtain 

The Poisson equation, also averaged over the horizontal plane and assuming RVl is independent of 
i?, is 

^=47rG((/,)+p.), (A3) 
where p-t is the background stellar density. 

If we now define a1 = (c^ + v1) and assume that this total velocity dispersion is independent 
of height z, we can combine the vertical momentum equation with the Poisson equation to obtain 
a second-order differential equation for the density profile (/?) p{z): 

d (-1 Mz),_ .4.G(p.+p(.)), (A4) 



dz \p{z) dz 

Henceforth, we assume that p^, is uniform within the gas disk, which is a good approximation 
provided that the stellar scale height is significantly larger than the gaseous scale height. Equation 
(IA4p is the expression for effective hydrostatic equilibrium in the vertical direction. 

Introducing a variable f{z) = ln{p{z) / p^) and a constant /i^ = (T^/(47rG/5*), we have 

f" = -^{l + ef), (A5) 
where the prime denotes a z derivative. This can be integrated once as 

^ = + + = (/o - / + - e^) , (A6) 

where /o = ln(po/p*) is the boundary condition at the midplane where /' = 0. If we Taylor expand 
and retain only the two lowest order terms, i.e. f{z) = fo — fiz'^ so that p/po = exp(— /iz^), the 
governing ODE becomes an algebraic equation: 

_ ( \ _ 47rG(po + p,) 2 f 
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The coefficient /i is 



Tlierefore, tlie gas density and pressure are approximately given by Gaussian profiles 



p{z)=poe 2H\ P{z)=Poe 2H\ (A9) 

where Pq = cr'^po- The midplane gas density po is determined by requiring that the profile integrates 
to the (known) gas surface density, 

p{z) dz = \/2tt pqH. (AlO) 



Substituting for pQ in equation (|A8p . the scale height H must satisfy 



2 

= (All) 



This yields a quadratic equation for H, with solution given by equation ^ of the text 



